Dynamical spin Hall conductivity in a magnetic disordered system 
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We investigate the intrinsic spin Hall effect in a quantum well semiconductor doped with magnetic 
impurities, as a means to manipulate the carriers' spin. Using a simple Hamiltonian with Rashba 
spin-orbit coupling and exchange interactions, we analytically compute the spin Hall conductivity. 
It is demonstrated that using the appropriate order of limits, one recovers the intrinsic universal 
value. Numerical computations on a tight-binding model, in the weak disorder regime, confirm that 
the spin Hall effect is preserved in the presence of magnetic impurities. The optical spin conductivity 
shows large sample to sample fluctuations in the low frequency region. As a consequence, for weak 
disorder, the static spin conductivity is found to follow a wide Gaussian distribution with its mean 
value near the intrinsic clean value. 

PACS numbers: 72.25.Dc, 72.25.Rb, 75.50.Pp 



I. INTRODUCTION 

The possibility to manipulate the electron spin by elec- 
tric currents was first proposed by Dyakonov and Perel, 
forty years agoii They noted that, in analogy with the 
anomalous Hall effect, an electric current must create a 
spin flow perpendicular to the charge current. This "spin 
Hall" effect is a consequence of the asymmetry, driven by 
the disorder potential through spin-orbit interaction, in 
the spin and momentum of the electrons scattered off 
impurities. Renewed interest in this effect arose more re- 
cently, when an intrinsic counterpart was predicted in 
strongly spin-orbit coupled semiconductors. ^'■i^i The in- 
trinsic effect results from the spin-orbit coupling in the 
underlying topology of the Bloch states (in momentum 
space), that gives rise to an anomalous component in 
the group velocity, perpendicular to the applied electric 
field. ''■^ It is based entirely on the one-particle band struc- 
ture, and therefore subject to the infiuence of disorder. In 
the case of an n-doped semiconductor quantum well, the 
intrinsic contribution to the spin Hall conductivity turns 
out to be exactly cancelled by the effect of scatteringi^""— 
The exact cancellation of the intrinsic spin Hall effect is 
particular to the Rashba spin-orbit coupling, linear in 
momentum, and is not present for instance in the Lut- 
tinger model of a p-type semiconductor.^^ 

Experimental observation of the spin Hall effect was 
achieved first by optical techniques to detect the spin 
accumulation,— and more recently by direct all- 
electrical measurements*^ The inverse intrinsic effect, 
that is the generation of a longitudinal charge current 
by the injection of polarized electrons has also been 
report edji^ 

The vanishing of the spin conductivity in the Rasbha 
model of a two-dimensional electron gas, for arbitrarily 
small disorder, can be traced back to the fact that the 
spin current operator = (h/A){a^,v}, is proportional 
to the spin s derivative 



where cr — {(Tx, Cy, az) is the vector of the Pauli matri- 
ces, V is the electron velocity, H is the Hamiltonian, and 
(x, y, z) refer to the directions of the applied electric field, 
the generated transverse spin current, and the out of 
plane spin polarization (Fig. [T]) . The commutator [H, Sy\ 
is not affected by the presence of randomly distributed 
impurities. In a stationary state, the time derivative of 
the averaged spin density vanishes, and hence the spin 
current must also vanish/^i^il? However, if the impurities 
were magnetic this argument is no longer valid, the spin 
current would depend on additional contributions, and 
could therefore remain finite. As a consequence, mag- 
netic impurities should in principle be able to restore the 
intrinsic Hall effect in a two-dimensional n-doped semi- 
conductor. 

The influence of magnetic impurities on the intrin- 
sic Hall effect was investigated using numerical compu- 
tations of the Landauer formula,^^ linear response the- 
ory and Kubo formula, and quasiclassical Green 
functions. It is well established that, contrary to the 
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FIG. 1. (Color online) Schematic geometry of the spin Hall 
effect; a longitudinal current driven by an external elec- 
tric field induces a transversal spin current that results in a 
polarization of the carrier spins following the sign of their 
momentum (spheres with arrows). 
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non-magnetic case, the spin Hall effect does not vanish 
in this case, although different values of the spin conduc- 
tivity in the clean limit were foundi^ Interest in mate- 
rials combining the electronic properties of semiconduc- 
tors with fcrromagnetism was enhanced by their poten- 
tial applications in spintronicsi^^— The main focus in 
the study of the so-called "diluted magnetic semiconduc- 
tors" is on the mechanisms governing the interactions 
between impurities, that determine their thermodynami- 
cal and magnetic properties)^ However, it would also be 
interesting to investigate the spin dependence of the car- 
riers response, to an applied electric field. Indeed, while 
the topological component of the anomalous Hall effect is 
widely documented in magnetic semiconductors, its spin 
current counterpart, the intrinsic spin Hall effect, is still 
a developing topic. This is the question we address in 
the present paper. Our main goal is to compute, using 
linear response theory, the intrinsic spin conductivity in 
a confined electron gas taking into account Rashba spin- 
orbit coupling and magnetic disorder. In particular, we 
want to obtain an analytical expression of the dynami- 
cal spin conductivity, to calculate first the clean limit for 
fixed frequency, and then the static limit. The order of 
limits is important here because the large scattering time 
limit does not commute with the small frequency limit »^ 
In Sec. ini we consider a simple model suitable to de- 
scribe the linear response of a confined electron gas to 
an external electric field, as determined by the spin-orbit 
coupling (inversion asymmetry) and scattering on mag- 
netic impurities (assumed to be in a paramagnetic state). 
We neglect extrinsic effects and non-magnetic impuri- 
ties, and take into account an effective mass conduction 
band. In Sec. IIIII we generalize the model to a finite 
band tight-binding Hamiltonian in a square lattice. As 
in the effective mass model, we consider spatially dis- 
tributed magnetic impurities with random orientations. 
We use a Chebyshev pseudo-spectral method to compute 
the density of states (total and local) and the spin Hall 
conductivity)^ Finally, we discuss and summarize our 
results (SecHVl. 



II. SPIN HALL DYNAMICAL CONDUCTIVITY 

In this section we compute the spin conductivity 
Cs_ff(w) as a function of the frequency w, starting with 
a simple model Hamiltonian H. The Kubo formula is 
analytically calculated in the weak disorder limit, taking 
into account the ladder diagrams. 



domly distributed magnetic impurities, 



H=^--rT-{z-Kp) + V{x) . 
Im n 



(1) 



where p and x are the momentum and position oper- 
ators of the electron in the plane {x,y). The impurity 
potential energy is assumed to be short ranged and char- 
acterized by an exchange interaction constant J, > and 
a microscopic length a. 



V{x) — JsO^ fii ■ cr 6{x — Xi) 



(2) 



where / is the set of impurity sites Xi, and fii the unit 
vector in the direction of the impurity magnetic mo- 
ment; we assume that the impurities are in a param- 
agnetic state, so fii = (sin^j cos </>!, sin^i sin^i, cos^i) is 
uniformly distributed over the sphere, with 9i G (0, tt) 
and € (0, 27r) independent random variables. Typi- 
cal orders of magnitude, for a III-V semiconductor het- 
erostructure, are: effective mass m = 0.067 me (me is the 
electron mass), lattice length a ~ 0.56 nm, Rashba con- 
stant A ^ 10^^^ eVm, and s-d exchange constant (for Mn 
impurities, for example) Jg 0.1 eV^^^ It is convenient 
to choose units such that h — m = a = e= 1. 

The clean Hamiltonian becomes diagonal in the chiral 
basci^ 
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where (fc, (p) are polar coordinates in momentum space 
kx + ifcj/ = . The corresponding eigenvalues are. 



e±(fc) = Y±Afc. 



(4) 



The transformation matrix between the spin and the chi- 
ral basis is 



(5) 



For a given Fermi energy ep, the spin splitting is given by 
A — 2Xkp, where kp — y/'^.tp is the Fermi wavenumber 
(energy is measured in units of fP/ma^). In the chiral 
base the velocity operator = i [H, x] is given by, 



{k — X)cos(p — iAsinyj 
iA sin if {k + X) cos (p 



(6) 



jx = —Vx is the current in the x direction, and the spin 
current operator j'^ = (1/4){ct2, Uj,} is, 



A. Model Hamiltonian 

We consider a two-dimensional system of noninteract- 
ing electrons of effective mass m, charge — e, and Rashba 
spin-orbit coupling constant A, in the presence of ran- 
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Note that the velocity operator, at variance with the spin 
current operator, explicitly depends on the spin-orbit 
coupling constant. We use the conventional definition 
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FIG. 2. (Color online) The conductivity is calculated in the 
non-crossing approximation, as a zero order loop with disorder 
averaged Green's functions, plus the ladder diagrams. The 
ladder diagrams are calculated by renormalization of the spin 
current vertex. 



of the spin current operator jy, but other choices are 
possible (in the presence of spin-orbit coupUng it is not 
a conserved quantity) 

In the presence of isotropic magnetic impurities, the 
scattering strength is characterized by a spin-flip relax- 
ation time 

- = 2ttNfcJ^ « cJf (8) 

T 

where c is the concentration of impurities (per unit area), 
and Np is the density of states at the Fermi surface; it 
approaches the value Np « l/27r when the Fermi energy 
is much larger than any other energy in the system {Np w 
TO/27rfi,^ in dimensional units) 



B. Kubo conductivity 

The spin Hall conductivity can be calculated in the 
framework of linear response theory;^ using the following 
Kubo formula,5S, 



U-U' {a\j'y\a'){a'\jx\a) 



OL.Ol' 



(9) 

where w is the frequency of the applied external elec- 
tric field, io is a small imaginary energy added to ensure 
causality, 1? is the system's area, to. and |a) are eigenval- 
ues and eigenstates of H, respectively, and /„ = /(cq.) is 
the Fermi energy function. Formula ^ is well suited for 
numerical computations of the tight-binding model^i^i 
(see Sec, mil below): for the conduction band model ([T|) it 
is convenient to express it in terms of retarded G^, and 
advanced G"^, Green functions,— 

aly{u) (Tr [j,^(fc)G«(e + uj, k)Uk)G<{e, k) 

+ j^{k)G<{e + u:,k)Mk)G^{e,k)]), (10) 

where 

_ , , f cPk f de ^ , 
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FIG. 3. (Color online) Real and imaginary parts of the intrin- 
sic spin Hall conductivity for a pure system p3[) -(|14 | ). The 
four branches ±aj± are indicated by vertical dashed lines, their 
separation is |a;+ — 4A^ (arrows). 



Tr^ is the trace over the spin index, (...) denotes the 
average over the spatial disorder and magnetic moment 
orientations, and G<(e) = /(e)(G^(e) - G^(e)). Calcu- 
lation of pUl) can be done using the diagrams of Fig. 
where solid lines represent the disorder averaged Green 
functions, (G^'^(e, A;))>ia 
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e-e-(fc) + 2^ 




V 



(11) 



and 



1 







V 



e-e+(fc) + ^/ 




1 



(12) 



both written in the chiral basis. The perturbation ex- 
pansion of Fig. [2] contains the lowest order one-loop con- 
tribution (the first diagram) and higher order one-loop 
contributions, containing disorder averaged terms, from 
which we only retain the non-crossing diagrams (ladder 
approximation) . 
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Before proceeding with the computation of these dia- 
grams, we present the spin Hall conductivity in the clean 
limit, for which the Kubo formula can be calculated ex- 
actly. The real part of cFsh is, 

4A2) - A2 



ReCTs//(w) = 



-1 
8^ 



1 - 



8A2 



log 



ujiuj - 4A2) - A2 



(13) 

which has branch cuts at the roots ±uj± with w± = 
2A2 ± V4A4 + A2 (Fig. O. This formula coincides with 
the result of Ref. |35| in the small frequency limit. The 
imaginary part is given by 

(CJ2_A2)2 



(4A2a;)2 



(14) 



where Q is the Heaviside function (see Fig. [31 bottom) . 

For the evaluation of cr = cr^y, using the diagram ex- 
pansion, 

of Fig. [21 we consider the weak scattering limit 1 /r, A <C 
ep, in which the products G^G^ or G'^G^, of like Green 
functions are small compared to the product G^G^: 

aly{u) = {Tr[l^^±^^l-l^jlG% + u.)j.G\e)])^ 

(15) 

where we omitted the momentum variable fc. Moreover, 
we consider the low frequency limit w <C ep, so the dif- 
ference of Fermi distributions is different from zero only 
in a small neighborhood of e^?. The calculation of the 
first digram is straightforward. 



-1 

8^ 

-1 



1 - 



(1 - irw) 



(l-i™)2 + A2r2_ 

t2 

r , W = 1 — ITLJ 



87rA2T2+u>2' ^ (16) 

where the first term gives the intrinsic contribution (t 3> 
1) to the spin Hall conductivity (— b/Stt in dimensional 
units). To this order, for which electron-impurity corre- 
lations are neglected, we obtain the same result as in the 
non-magnetic case.-^^ 

Next, we compute the vertex correction (ladder dia- 
grams, and third line in Fig. [2]) in the same approxima- 
tion, where the Fermi energy e^^ , is much larger than any 
other energy scale: 

a^H{^) = (YT[ry{u)G"ieF+oj)j^G^{eF)]) (17) 

where the trace now runs over momentum and spin, and 
Jy = Jy{uj, k) is the corrected spin current vertex, which 
is a function of the energy uj). The algebraic equation 
that determines the corrected vertex, the third line of 
Fig. [21 is conveniently expressed in the spin basis, 

kdkdip f sm6d6d<l) f cosd e~'''^sin( 
I e'"^ sin 9 - cos 6* 



V — . 



(27r)2 J An 
X Uiip)G^iO,k)UH^) [j^' 
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(18) 
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FIG. 4. (Color online) Dynamical conductivity asH as a func- 
tion of UJ, for T varying between 40 and 190 (light to dark 
lines). (Inset) As — > the real part tends to l/8n for 
UJ > 1/r, but for cj < 1/t it increases up to 8/7 « 1.14. 
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FIG. 5. (Color online) Static spin Hall conductivity as a func- 
tion of the disorder strength from Eq. (|22[) . for two values of 
A. The dot at = corresponds to the clean static limit. 
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where we used the unitary matrix (O; we note that Jy, 
in the spin basis, has the structure, 



J" 



with J^_i = J;^ given by, 



J- 



^ti 



ikpATw 



n 



2«;2 + 2(7w + irw)(A2r2 + ^2) 



(19) 



(20) 



Introducing this expression into (|17p and computing the 
trace, resuhs in the ladder contribution. 



- AlTUJ 



+ (7W + iTtj)(A2T2 + w2) 



(21) 



Cohecting the two contributions (fTOI and (pij) we fi- 
nally obtain the dynamical spin Hall conductivity in the 
presence of magnetic impurities, in the weak scattering 
strength approximation, 



CTsH (i^) 



-1 



StT A2t2 



1- 



AlTUJ 



2u;2(4 - 3i™) + A2t2(7 - Giro;) 



(22) 



This result is similar to the one obtained in Ref. 21, al- 
though it has an extra term 4irw in the numerator of the 
last term; this difference, that does not change the limit- 
ing behavior for small frequency or large scattering times, 
can be attributed to the different approximation methods 
used. We show in Fig. [¥]the real and imaginary parts of 
UsH for different values of the disorder scattering time, 
as a function of the frequency. The real part of UsH, an 
even function of the frequency, has a sharp change of sign 
at w = A in the clean limit, that smoothes out with finite 
T. The imaginary part shows, at the same frequency, a 
peak that broadens with increasing disorder. Comparing 
with Fig. [21 we note that the in the weak scattering limit 
the splitting w- + a;+ is neglected. 

The static spin current, as deduced from the real part 
of (1121) at w 0, 



fT.ff(O) = 



-1 8A2 



7A2r2 ' 



(23) 



is shown in Fig. O The continuous band model predicts 
an enhancement of the spin conductivity near the clean 
limit T ^ A^^, and a slow reduction for smaller scatter- 
ing times. 



C. The clean limit of asH 

The imaginary part of Ush is an odd function of the 
frequency, and thus it must vanish at zero. Moreover, 
the real part, contrary to the non-magnetic case, does 
not vanish in the clean limit. For t — > oo wc find. 



lim aly{uj) = 



-1 A2 
Svr A2 - 



(24) 




FIG. 6. (Color online) Convergence of the conductivity for 
values of uj near zero, ranging from 10~^ to 10~^. The ex- 
pression of 8irRe [o"(t<-')] for u = 10~" converges to zero only 
for r > 10"^^. The horizontal line in the graph corresponds 
to the numerical factor 8/7. 



if now we put w — >■ we recover the intrinsic value of 
the static spin Hall effect, as might be expectedjSi How- 
ever, the dynamical spin conductivity near the origin is 
singular, 29 as can be inferred from the behavior of Re asH 
near oj ~ (inset in Fig.jH). Indeed, we observe that con- 
vergence to the static limit is not uniform: the larger the 
value of T, the steeper the slope at the origin. This is 
related to the fact that the limits, r — > oo and a; — do 
not commute: 



lim aUO) = 



-18 

' 



(25) 



as seen from formula ([231). Various different results for 
the clean limit of <Jsh (0) were reportedj ^^i^9'2^ The reason 
can be attributed to the fact that the static value cTsh (0) 
cannot be computed at a; = 0, for r tending to infinity. 
The order of limits starting with a; — > 0, would violate 
the causality condition contained in the linear response 
formula. 

We present in Fig. [51 a sequence showing the spin con- 
ductivity as a function of the scattering time for the se- 
ries ujn = 10^", with n = 2, . . . , 9. We observe that, for 
a fixed value of the frequency a;„, convergence to — I/Btt 
is reached for r > r„+i (where r„ = 10"). Therefore, 
the limit w — > does not converge for any finite value of 
T. In fact, one can erroneously interpret the transitory 
value of 8/7 as being the static spin conductivity. 

We conclude that in the presence of magnetic impuri- 
ties, and in the weak disorder approximation, the static 
spin conductivity jumps from its intrinsic value to a value 
about 8/7 for finite disorder, and then it smoothly de- 
creases when the disorder increases. 
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III. THE TIGHT BINDING MODEL AND 
NUMERICAL RESULTS 

We will now go beyond the effective mass model and 
the perturbation approximation in the disorder strength, 
to compute the spin conductivity. We introduce a dis- 
crete tight-binding model to investigate, in particular, 
the influence of a finite band width on the spin transport. 
The Hamiltonian is defined on a finite square lattice of 
size L X L and lattice constant a, with Rashba spin-orbit 
coupling A, and exchange interaction on impurities Jg, 



H = -tY^ c\cj + - ^(i4cr:rc,+.g 
-I- Js'^c^^hi ■ <jci, 



(26) 



where the hermitian conjugate terms of the first two sums 
are implicit; c\ — (cj^ c\^) is the creation operator at site 
I, of spin up t or down J,, and is the corresponding an- 
nihilation operator. We denoted x and y the two orthogo- 
nal directions in the lattice. The kinetic term is limited to 
nearest neighbours and its bandwidth is 4t. This model 
is a straightforward generalization of ([1]), with discrete 
kinetic, spin-orbit and exchange terms. In the following 
we take i = 1 as the energy unit, and fi = a = e= las 
before. 

The computation of the Kubo formula (O needs in 
principle the knowledge of the whole spectrum and the 
full set of eigenvectors of H . One can avoid this numeri- 
cally expensive task, by reducing the formula to a trace, 
which can be evaluated with a stochastic method, in any 
basis. Specifically, we adopted a pseudospectral imple- 
mentation of the kernel polynomial method.'*^ 

Because the real and imaginary parts of the spin con- 
ductivity ^ are related by the Kramers-Kronig rela- 
tions, it is enough to compute the optical spin conduc- 
tivity (To = Im UsH , defined by, 

^o(^) = E a')^^^^S{u; - e„ + e„0 , (27) 

a,Q' 

where 

J^y{a,a')^lm[{a\fy\a'){a'\j^\a)] , 
can be written in the form, 

aoiio) ^ - fde j{e + e) [/(e) - f{e + u)] , (28) 

using the definition, 

a, a' 

The advantage of this representation is that the matrix 
j(e, e') is a function of energies only which can be ex- 
panded in Chebyshev polynomials products T„(e)Tm(e') 



with expansion moments /x„m, which are easily computed 
using a stochastic trace evaluation method. We now 
briefiy detail the algorithm. 

The first step is to scale and translate the energy, in 
such a way that the Hamiltonian spectrum is mapped 
into the interval e G [—1, 1]. The series expansion of 
can then be written. 



n,m— 
M-1 

E 

n,m— 



(30) 



where Af^ is the number of moments we use, T,„ is the 
Chebyshev polynomial of order m, 

9 9 

hnm. — 



1 + Sn,o 1 + Sjn,0 ' 



and gm is a filter that minimizes the Gibbs oscillations 
arising in truncating the series to a finite order (we use 
the Jackson kernel, see Ref. IsgI ) . 

The second step is to compute ^inm- The presence 
of the radicals in the denominator allows to express the 
moments as usual scalar products. 



^J■r. 



1 /.I 



dede'j{e,e')T,,{e)T,n{e') (31) 



-1 J -1 

which after rearranging, 

Mnm = ■^lm^{a\jy\a'){a'\ja:\a)Tn{ea)Tjn{e'a') 



a, a' 



^lmJ2{a\T^{H)j;\a'){a'\T^{H)j,\a) 



a, a' 



can be written as a trace formula, 
1 



ImTr[r„(i?)j-r„,(iJ)j, 



(32) 



In the third step, the trace is evaluated using an average 
over random states |r), 



-E 



;r|T„(i/)j-T„(iJ)j,|r), 



where N^. is the number of random vectors needed for av- 
eraging. Next, one uses the recurrence relations between 
the Chebyshev polynomials to obtain the operator prod- 
ucts to order (n, m). In our case, even if the spin current 
operator is traceless, fast decrease of statistical errors 
and a rapid convergence of the algorithm is obtained. In 
order to optimize the Chebyshev recursion algorithm we 
adapted a pseudospectral procedure: 

H\r) = IFFT[i/o(fc)FFT[|r)]] + V\r) , 
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where Hq contains the kinetic and spin-orbit terms, di- 
agonal in the /c-representation, V is the impurity poten- 
tial, diagonal in the x-representation, and FFT, IFFT 
denote direct and inverse fast Fourier transforms, respec- 
tively. Once fimn are known, we can efficiently obtain 
j(e, e') from formula using collocation points 
e; = cos[7r(/ -f- l/2)/Me]. One drawback of this method 
is that the cj-dependent integral (|28p . cannot be directly 
computed, since e; -I- w is not necessarily a collocation 
point. One therefore has to perform an interpolation of 
some sort. 

The last step is thus to compute the convolution inte- 
gral (pS)) . We insert into this equation a (5-function in its 
integral form, in order to separate the energy integrals. 
The expression for (Tq becomes, 



(To(w) 



dt 
2^ 



dede 



,[£ /leit(.-.'+-)j(£^g') 



We discretize the frequencies as well as the energies and 
the auxiliary time variable t. The sampled ao, finally 
takes the form of a Fourier transform 



1 



where 



and so ao{oJn) is straightforwardly evaluated by a fast 
Fourier transform of F . From the expression of the op- 
tical spin conductivity we calculate its static value using 
the Kramers-Kronig relation. 



The density of states. 



dio o'o(w) 



(33) 



(34) 



may be computed using the same method, and the ex- 
pansion 



M-l 



Mm hfnffmTn (^) 



where 



IJ-r, 



In order to avoid finite size artifacts and to approach 
the thermodynamic limit faster, we use twisted boundary 
conditions: ^{x + L,y) = e*'^^5'(a;, y) and "^{x, y + L) — 
e*¥'H\[/(x, y), for any state Physical quantities are aver- 
aged over many phase configurations {(px,(py) G [0, 27r]^. 
This method has already been applied to determine the 
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FIG. 7. (Color online) Density of states for cj^ — (clean 
system), and cjf = 0.01, 0.06; A = 0.4. The spin-orbit in- 
teraction is responsible for the splitting of the chiral states 
(central peaks). 



finite-size spin Hall conductivity, which in the case of 
non magnetic impurities, should vanish in the thermody- 
namic limit 1^1^ In our case, the spin Hall conductivity 
being finite even at the thermodynamic limit, twisted 
boundary conditions constitute a useful tool to improve 
numerical convergence. For instance, the translational 
invariance of the disorder in the limit of infinite number 
of configurations, is recovered by averaging over twisted 
boundary conditions. 

Let us make some remarks about the choice of the num- 
ber of moments M. Filtering the truncated Chebyshev 
expansion, introduces a broadening of sharp features that 
scales as causing a reduction in the energy resolu- 

tion. The optimal number of moments can be calibrated 
from the clean (Js — 0) case, and it turns out to be a 
few times the linear system size M = 2L, 4L, therefore 
4L^, 16L^ coefficients are involved in the j(e, e') expan- 
sion. We obtained a a well detailed density of states, 
without finite-size artifacts. This choice of numerical pa- 
rameters, has been verified by comparing small samples 
averaged over twisted boundary conditions, with equiva- 
lent samples of greater size. 



A. Numerical results 

In the following, we describe the properties of spin 
transport through the numerical computation of the den- 
sity of states, and the static and optical spin conduc- 
tivities. Simulations were performed on lattices of size 
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FIG. 8. (Color online) Optical spin conductivity as a func- 
tion of the frequency, for different disorder strengths. For 
each = cjg — 0.002,0.01,0.06, three values of the con- 
centration c = 0.04, 0.06, 0.08, were used; eF = 2, A = 0.2. 



64^, 128^, but higher sizes were used to test convergence. 
Typical numerical parameters are: A = 0.2, 0.4, in units 
of at, c from 2% to 10%, and t from 5 to 500, in units of 

h/t. 

We will limit the range of parameters in order to pre- 
serve the qualitative shape of the clean density of states, 
avoiding for instance an impurity band regime (strong 
disorder regime). This choice should allow us to compare 
the tight-binding model to the analytical results of the 
preceding section, and to determine its validity. One ob- 
serves in Fig. [7]that, within the range = 0, . . . , 0.06, 
the density of states p{e) remains close to its clean limit, 
and the spin-orbit splitting, proportional to A, fades out 
with stronger scattering (to enhance the effect of the spin- 
orbit coupling, we put A — 0.4). 

The spin conductivity depends on the Fermi energy 
ep, the spin-orbit coupling A, and in principle, the con- 
centration c and the exchange constant Jg. In fact, the 
scattering off impurities is characterized by the sole pa- 
rameter T"^ = cJ^, at least in the weak disorder limit 
<C 1, hence we would expect that (Tq — o'o(w;r). In 
Fig. [5] we represented the optical spin conductivity as a 
function of the frequency, for different sets of pairs c and 
Js, giving the same values of r. We chose the Fermi en- 
ergy at the middle of the band, ep = 2. We observe that 
the curves coincide according to their t values, confirm- 
ing that in the regime studied, the relevant parameter 
is the product cJ^ . The maximum of ao corresponds to 
w « A as in the continuous case, and is almost indepen- 
dent of T (compare with Fig. [3]). In Fig. [21 we present (Tq 
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FIG. 9. (Color online) Optical spin conductivity as a function 
of frequency and ef, for two different disorder strengths = 
0.002, 0.2. As the disorder increases, the optical conductivity 
decreases and spreads out over uj. 



as a function of the frequency and the Fermi energy. Its 
maximum appears at Fermi energies in the middle of the 
band, and follows a law in w ^ v^cf, compatible with the 
estimation w « A given above. The intensity of the peak 
decreased by a factor 2, for a factor of 100 in t (it is 500 
for the top figure and 5 for the bottom one). The optical 
spin conductivity vanishes for small a; <C A, and large 
frequencies w ^ A. Therefore, the effect of increasing 
the disorder is, as expected, to broaden and damp the 
peak on the spectrum. However, at variance with the 
conduction band model for which the peak width is pro- 
portional to T~^, in the tight-binding model, the effect 
of a finite band width quantitatively change the scaling 
with the scattering time (the large Fermi energy limit is 
non longer applicable). 

The computation of the static spin conductivity (j33p . 
appears to be particularly sensitive to statistical errors, 
as already noted^Si A large set of boundary conditions 
and random impurities distributions, is necessary. Re- 
sults are shown in Fig. IIOI where a^H (0) is represented as 
a function of the Fermi energy and the disorder strength. 
The spin Hall conductivity is an odd function, it vanishes 
at ep = and outside the energy band, but is nearly con- 
stant over a large range of Fermi energies. Its slow fall 
off with disorder underlines the weakness of vertex cor- 
rections in the magnetic case, as shown by the analytical 
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1.0 




FIG. 10. (Color online) Static spin Hall conductivity as a 
function of ep, for disorder strengths — 0, 0.2. (Inset) 
Integrand of H33|) . used to compute (7sh(0); the error bars 
represent the amplitude of cro(i^) /i^ fluctuations near the zero 
frequency. 

calculation of the previous sectioni^ More importantly, 
we observe large fluctuations of the spin conductivity at 
low frequencies, depending on the sample distribution. 
Strong spin conductance fluctuations are known to arise 
at mesoscopic scalesjS and, provided that scattering en- 
sures ergodicity, their universal statistics can be related 
to that of random matricesi^ We found a noticeable am- 
plification of the fluctuations in the region uj < that 
explains the large dispersion of asH (0) values around its 
mean. In the inset of Fig. [SJ the integrand CTo/o; in 
is represented together with error bars proportional to 
the sample to sample dispersion, 

Aa,^(o.) = {a,H{io)' - {a,H{io))Y'^ , (35) 

where the angle brackets stand for the random potential 
averaging. The averaged static spin Hall conductivity, 
is then not enough to characterize the spin transport. 
In fact, even in the present weak disorder limit, where 
localization of the carriers can be neglected, we find a 
Gaussian distribution of the spin conductivity fluctua- 
tions, as illustrated in Fig. [TlJ The standard deviation 
of the static spin conductivity Acrs-ff(0)/(crs_y (0)) ~ 0.1, 
giving an incertitude of 10% about the mean value. 

IV. CONCLUSION 

In this paper we investigated the influence of magnetic 
impurities in the intrinsic spin Hall effect. A first model. 
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FIG. 11. (Color online) Normalized histogram of the static 
spin conductivity for = 0.2, showing a Gaussian distribu- 
tion with mean value {(Tsh{'^)) = 0.98 [— I/Stt] and standard 
deviation Acrs£r(0) =0.1 (solid line). 



amenable to analytical calculations, was introduced. We 
considered a simple one conduction band Hamiltonian in- 
cluding the effects of Rashba spin-orbit coupling and ex- 
change interaction with randomly distributed magnetic 
moments. Using standard linear response methods we 
obtained, in the weak scattering and low frequency ap- 
proximation, the expression of the dynamical spin con- 
ductivity as a function of the system's parameters (spin- 
orbit coupling, scattering time, Fermi energy). We em- 
phasized the relevance of the order of the large scat- 
tering time and small frequency limits, and found that 
when frequency is going to zero, one must recover the 
spin conductivity of an ideal two dimensional electron 
gaz, a — — e/87r. However, an arbitrarily small amount 
of magnetic impurities will increase the static spin con- 
ductivity by a factor of about 8/7, showing the singular 
nature of the clean limit. 

We also investigated, using numerical methods, the 
tight-binding generalization of the conduction band 
Hamiltonian. We adapted the kernel polynomial method 
using a pseudospectral evaluation of the Chebyshev re- 
cursion to compute the density of states and the opti- 
cal spin conductivity. As a function of frequency, the 
spin Hall conductivity shows a peak at a; ~ X^/tF, that 
broadens and weakens with the increase of the disorder 
strength. These changes are determined by the combina- 
tion ~ cJg , that characterize the disorder strength. 

Within the regime of weak disorder, that is in the ab- 
sence of separate impurity bands, the deviation with re- 
spect to the clean system is small. The main effect of 
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magnetic impurities in contrast with non-magnetic ones, 
is to preserve and reinforce the intrinsic mechanism that 
would otherwise be destroyed. In the finite band model, 
we observe an amphfication of spin conductivity fluctu- 
ations at low frequencies, which are at the origin of a 
significant incertitude of the static UsH mean value. A 
detailed theory of the spin Hall conductivity in the pres- 
ence of magnetic impurities should consider these sam- 
ple to sample fluctuations, that could become the main 
physical effect in the description of spin transport in the 
strong disorder regime. 

In spite of the absence of magnetic order (we treated 
the paramagnetic state), the sole presence of spin-flip 
processes, guarantees the emergence of a spin current 
proportional to the applied electric field. For weak disor- 
der, and strong spin-orbit coupling, the intrinsic spin Hall 



effect is reinforced by multiple scattering interactions, de- 
scribed by the vertex correction to the spin current. Ma- 
nipulation of the carrier's spin in a semiconductor impu- 
rities, is then possible by doping it with magnetic impu- 
rities and applying an electric field. This result naturally 
raises the question of the interplay between intrinsic spin 
Hall currents and ferromagnetic impurities as present in 
a diluted magnetic semiconductor. 
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